Restricting UHECRs and cosmogenic neutrinos with Fermi-LAT 
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Ultrahigh energy cosmic ray (UHECR) protons interacting with the cosmic microwave background 
(CMB) produce UHE electrons and gamma-rays that in turn initiate electromagnetic cascades on 
CMB and infrared photons. As a result, a background of diffuse isotropic gamma radiation is 
accumulated in the energy range E < 100 GeV. The Fermi-LAT collaboration has recently reported a 
measurement of the extragalactic diffuse background finding it less intense and softer than previously 
measured by EGRET. We show that this new result constrains UHECR models and the flux of 
cosmogenic neutrinos. In particular, it excludes models with cosmogenic neutrino fluxes detectable 
by existing neutrino experiments, while next-generation detectors as e.g. JEM-EUSO can observe 
neutrinos only for extreme parameters. 



I. INTRODUCTION 

The origin of ultrahigh energy cosmic rays (UHECRs) 
is not yet established despite more than 50 years of re- 
search. Natural candidates as UHECR primaries are ex- 
tragalactic protons from astrophysical sources. In this 
case, interactions of UHE protons with the cosmic mi- 
crowave background (CMB) leave their imprint on the 
UHECR energy spectrum in the form of the Greisen- 
Zatsepin-Kuzmin (GZK) cutoff and a pair-production 
dip Q. 

The GZK cutoff is a steepening of the proton spectrum 
at the energy Eqzk ~ (4 — 5) x 10 19 eV, caused by photo- 
pion production on the CMB. Such a steepening has been 
observed by the HiRes @ and the Auger collaboration Q , 
but its real cause is still unclear. 

An immediate consequence of the dominance of extra- 
galactic protons in the CR flux and their interaction with 
CMB photons is the existence of ultrahigh energy ( "cos- 
mogenic") neutrinos produced by charged pion decays, 
as suggested first in Ref. [4|. 

Another signature for extragalactic protons is a pair- 
production dip [fl Q in the CR flux around 5 x 10 18 eV, 
which is clearly seen in the experimental data. Photons 
and positrons from pion decay and p+jcmb p+e + +e~ 
pair-production initiate electromagnetic cascades on pho- 
tons from the CMB and the extragalactic background 
light (EBL), dumping all the energy injected into cas- 
cade particles below the pair-production threshold at 
-100 GeV. 

Clearly, the production of neutrinos by UHE protons is 
thus intimately tied to the one of photons and electrons, 
and both depend in turn on the flux of primary cosmic 
rays. While UHE photons and electrons start electromag- 
netic cascades by scattering on photons from the EBL, 
neutrinos reach us suffering no collisions. Therefore, 
the measurement of the diffuse extragalactic gamma-ray 
background (EGRB) can be used to impose a strict upper 
limit on the possible diffuse high energy neutrino flux, as 
suggested first in Ref. Q. 



We derive in this work an upper limit on the flux of cos- 
mogenic neutrinos assuming that the primary UHECR 
particles are protons. In case that all primaries or part 
of them are nuclei, the cosmogenic neutrino flux is lower 
than for a pure proton composition Q. Thus our as- 
sumption of a pure proton composition is justified, since 
we aim at deriving an upper limit on the cosmogenic neu- 
trino flux. Note that the HiRes data [2J agree with a pure 
proton composition at E > 1 X 10 18 eV, while the mass 
composition deduced from the Auger data indicates the 
presence of heavier nuclei in the primary UHECR flux 0] . 
In the latter case, the maximally allowed cosmogenic neu- 
trino flux would be below the upper limit derived for a 
pure proton primary flux in this paper. 

In the present work, we use a recently reported mea- 
surement @ of the EGRB by Fermi-LAT to constrain 
UHECR models. We show that the observed fast de- 
crease of the EGRB with energy, J(E) oc E~ 2A1 , already 
constrains such models. In particular, versions of the dip 
model with strong redshift evolution contradict the Fermi 
data, while this model without or with weak redshift evo- 
lution remains viable. Moreover, the Fermi data allows 
us to derive a strong upper limit on the diffuse UHE neu- 
trino flux. As a result, we conclude that the detection of 
cosmogenic neutrinos requires to increase the sensitivity 
of UHE neutrino experiments compared to current lev- 
els. As it is demonstrated below, the maximal energy 
density of cascade radiation w™f x ~ 5.8 x 10 -7 eV/cm 3 
allowed by the Fermi-LAT data can be used to select vi- 
able UHECR models without explicitly calculating elec- 
tromagnetic cascade processes. 



II. ANALYTICAL CALCULATIONS 

The two basic processes driving an electromagnetic 
cascade are pair production (PP) 77b — > e + e~ and in- 
verse Compton (IC) scattering e ±/ yb — > e ± 7 on back- 
ground photons 7b- The cascade develops very fast with 
a minimal interaction length h n t{E) — lOkpc until it 
reaches the pair creation threshold. From that point on, 
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electrons emit photons in the Thomson regime while pho- 
tons stop interacting. Their spectrum can be estimated 
analytically 0, E3] in terms of the production rate of the 
cascade photons Q'i a,s (E) per unit volume as 



K{E/e x )-^ 
K(E/e x )- 2 



for E < ex , 

for e x < E < e a . 



(1) 

with a steepening at E > e a . Here, e a is the minimal 
absorption energy of a cascade photon scattering on the 
EBL, and ex is the energy of a photon emitted by an 
electron/positron [ex + 7 — > el + 7.x"), which is in turn 
produced by a photon 7 a (via 7 a -I- 7ebl — > e x + e x ) 
with the minimal absorption energy e a . The energy spec- 
trum (fTJ) of the cascade radiation typically extends up 
to ~ 100 GeV. The constant K in Eq. © defines the 
normalization of the production rate via K = Q^ s (ex)- 
The two energies e a and ex are related to each other as 
£x = 1/3 (e Q /m e ) 2 e cmb 0,113], where e cmb = 6.35 x 10~ 4 
eV is the mean energy of CMB photons. 

We can account for the absorption of cascade photons 
on the EBL, integrating their production rate Q^ S (E) 
over the volume of the universe, 



)cas 



■ exp 



knt(E) 



(2) 



Integrating the rate from r = up to cH^ 1 we obtain 
the relation between the absorbed flux J^(E) and the 
unabsorbed flux J^ as (E), 



J^ S (E) = J™(E) 



knt(E) 



1 — exp 



cm 



knt(E) 



(3) 



Here, Ho denotes the present value of the Hubble param- 
eter. 

The cascade energy density u; cas at the present epoch 
is calculated as 



^ JdEEJtlUE) 



(4) 



In Fig. [TJ we show the measurement of the EGRB by 
Fermi-LAT Q (black circles with error bars) together 
with the maximally allowed photon flux (solid red line) 
derived analytically. More precisely, we have determined 
the maximally allowed photon flux requiring that the 
curve just touches the lower end of the error bars of the 
Fermi-LAT data. The corresponding bound on the cas- 
cade energy density is w™f x = 5.8x 10~ 7 eV/cm 3 . 

The bound on w™| x derived by us can be used to select 
in a simple way viable UHECR models. 

We limit our consideration to pure proton-composition 
models, which are described by the generation index a g , 
the maximum acceleration energy E max and the cosmo- 
logical evolution of the sources parametrized by (1 + z) m 
with fixed m and maximal redshift z max . The quantity 
of interest, the space density of protons n p {E, z) at each 
cosmological epoch, is calculated as in Ref. [6| in the con- 
tinuous energy-loss approximation. To evaluate the role 



of fluctuations in» + 7 — > 7r+all scattering, the density of 
protons in Ref. [6| was computed also solving the kinetic 
equation (16). It was found that for E max = 1 x 10 21 eV 
both methods give practically identical results, while for 
E m &x = 1 x 10 23 eV the difference does not exceed 15%. 
In the calculations of this work, the UHE proton fluxes 
at z — are normalized to fit the HiRes spectra 0]. 

The UHE diffuse neutrino flux at highest energies de- 
pends mainly on ct g and -E max . The generation index is 
limited as 2.0 < a g < 2.7, by the following reason: Let us 
choose first the minimum possible index a g — 2.0. In this 
case, the calculated extragalactic UHECR flux is very 
flat and can explain the observed spectrum only above 
E ~ (0.5—1) x 10 19 eV, i.e. above the ankle. Increasing a g 
decreases the predicted transition energy between Galac- 
tic and extragalactic UHECRs, until for a g w 2.6 — 2.7 
this energy becomes lower than 1 x 10 18 eV, where, as 
observations show, heavy nuclei dominate. 

UHECR models allowed by Fermi-LAT data, i.e. lead- 
ing to a cascade energy density u; cas below u;™| x = 
5.8 x 10 _7 eV/cm 3 , are characterized by a low value of 
E ma ,x or weak source evolution. For each given model, 
we calculate w cas as 

w cas = / ^^-E(3 0tCm [(l + z)E}n p (E,z), (5) 
J 1 + z 

where n p is the (physical) density of protons at redshift 
z, Po(E) = (l/E)(dE/dt) is the relative rate of energy 
loss of a proton with energy E at z = 0, and /3o,em de- 
notes the relative rate of energy injected by protons into 
electromagnetic cascades due to pair production and pion 
production (pj — > 7r ± — > and P7 — > ir° — > 7) at z = 0. 
In Table 1, we report the numerical values obtained for 
w cas for different values of the maximal energy of accel- 
eration -E max , exponents a g of the generation spectrum 
and the maximal redshift z max . The UHE proton fluxes 
at z = are normalized to the HiRes data 0]. In Ta- 
ble 1 we show also the ratio of the contribution from 
pair production to the one from pion production. The 
cases with different a g in Table 1 correspond to different 
transition energies from galactic to extragalactic cosmic 
rays, increasing from E ~ 1 x 10 18 eV for a g = 2.6 to 
£~5x 10 18 eV for a g — 2.0. Table 1 gives examples of 
UHECR models that are allowed and that are forbidden 
by the Fermi-LAT data. 

The upper part of Table 1 (no-evolution case) presents 
the allowed models with w cas < 5.8 x 10~ 7 eV/cm 3 . In 
the two lower parts the cosmological evolution of UHECR 
sources is included, assuming that the product of the co- 
moving source density n s and the source luminosity L s 
evolves as n s (z)L s (z) = noLo(l + z ) m ■ The middle part 
contains both allowed and marginally allowed evolution- 
ary models. The large neutrino fluxes at highest ener- 
gies, favorable for detection by the JEM-EUSO instru- 
ment and by radio methods, are expected in the models 
with flat spectra (a g = 2.0) and large fi max . Examples 
of such models allowed by w cas are also given in Table 1, 
most notably the one in the first row of the middle part. 
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TABLE I: The energy density oj cas of the cascade radiation 
produced by UHE protons normalized to HiRes data. 



No evolution, allowed models 


m 


Emax 


«s 


Zm ax 


t^cas [eV/cm 3 ] 


e~*~ e i 7r 
^cas / ^cas 





10 21 


2.0 


2 


4.0 • 10" 8 


2.41 





10 21 


2.7 


2 


1.2 • 10~ 7 


24.5 





10 22 


2.0 


2 


5.1 • 10" 8 


1.2 





10 22 


2.7 


2 


1.2 • 10~ 7 


20.8 





10 22 


2.0 


3 


5.5 ■ 10" 8 


1.2 





10 22 


2.7 


3 


1.4 • 10~ 7 


22.0 


With evolution, allowed models 


2.5 


10 22 


2.0 


4 


2.7- 10~ 7 


1.46 


3 


10 21 


2.5 


2 


4.3 • 10~ 7 


14.7 


3 


10 22 


2.4 


3 


5.8- 10~ 7 


8.4 


Models excluded by ui caa 


3.5 


10 23 


2.3 


3 


7.4 ■ 10~ 7 


4.9 


4 


10 21 


2.5 


2 


8.6- 10~ 7 


15.3 


4 


10 22 


2.0 


3 


7.6- 10~ 7 


1.5 



On the other hand, large neutrino fluxes at energies up to 
1 x 10 17 eV, which are favorable for IceCube detection, 
do not require large E max but strong evolution. One 
can find such models in Table 1, too, among the allowed 
models. In the lower part of Table 1 three examples for 
evolutionary models forbidden by the Fermi-LAT data 
are shown. 

Note that apart from the cascade radiation, one should 
expect various additional contributions to the measured 
Fermi flux and these contributions lower w™f x fur- 
ther. Among them are photons from unresolved fill ] or 
dead [12| active galactic nuclei (AGN) and from other 
galaxies. Dark matter annihilations and/or decays in the 
extended DM Galactic halo or beyond can give another 
contribution to the Fermi flux. Subtracting these pro- 
cesses would strengthen the upper limit on the neutrino 
flux derived below, and thus our result is conservative. 
On the other hand, extragalactic magnetic fields with 
strength above 1 nG play the opposite role: High-energy 
cascade electrons loose energy radiating predominantly 
synchrotron photons, which are not detected by Fermi- 
LAT and thus the energies of these electrons do not con- 
tribute to w cas . We show below that in the case of mag- 
netic field strengths less than 1 nG this correction on the 
limit for UHE neutrino flux is small. 



III. MONTE CARLO SIMULATION 

In addition to the analytical treatment, we obtain the 
EGRB spectrum based on a Monte Carlo simulation of 
the cascade development. We generate CR sources from 
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FIG. 1: Fermi-LAT EGRB spectrum (black circles with with 
error bars) in comparison with maximally allowed fluxes given 
by analytical (solid red line) and MC calculations (red stars 
for ot g = 2.6, and blue stars a 9 = 2.0). All three curves are 
normalized by the highest energy point of the Fermi spectrum. 
The MC EGRB fluxes are calculated for the following values 
of the parameters: £ max = 10 21 eV, z max = 2, m = and 
Oi g — 2.6 or 2.0. Also shown two MC spectra with the same 
parameters as above, but normalized to the HiRes proton 
spectrum (red boxes for ct g = 2.6 and blue boxes for a g — 2.0 
) . The plot illustrates the universality of the cascade spectrum 
and reasonably good agreement between MC and analytical 
results. 



a homogeneous source distribution up to a maximal red- 
shift z max . Assuming the proton injection spectrum in 
the form dN/dE oc E~ ag, d(E — E max ), we propagate the 
UHE protons accelerated in the sources through the ex- 
tragalactic space, using the Monte Carlo code described 
in [13j], until their energy is below the threshold for e + e~ 
pair production, E m - m m 10 18 eV, or until they reach the 
Earth. For the simulation of pion production we use 
SOPHIA [T3|, while e + e~ pairs are injected according 
to the continuous e nerg y losses and their mean energy 
calculated in Ref. @, [l5[ 

We follow the evolution of electromagnetic cascades 
using the Monte Carlo code introduced in Ref. [l6| and 
the best-fit model of for the EBL energy density. 
The MC procedure provides an one-dimensional descrip- 
tion of the cascade development, taking into account the 
pair production and IC processes as well as adiabatic en- 
ergy losses. Extragalactic magnetic fields with average 
strengths close to the upper limit B ~ 1 nG have a small 
influence (of order 20%) on the resulting EGRB. We will 
discuss further this result below. 

In Fig. Q] the cascade fluxes are shown for two UHECR 
models. The curve marked as a g — 2.6 (red boxes) gives 
the cascade flux for the non-evolutionary (m = 0) dip 
model Q with E max = 1 x 10 21 eV and z max = 2 normal- 
ized to HiRes data. The other curve marked as a g = 2.0 



4 



10' 



-icr 



w 
\n 

E 

o 

> 

0) 

LU 10° 
LU 



10" 



: 


z =2, E =10 21 eV 






max max 


HiRes 


: 


*+l Fermi LAT 


V 


: 

r 


+ % 

T 




■ 


B 


P\ 


r y .- 






















/ T : 






/ v°°= " 






/ f; : 


j 


.J M..J M..J M..J M..J M..J /....J ,,„ji 


'■■J ...J J ....J .AJ 



10 6 10 8 10 10 10 12 10" 10 16 10 18 10 20 10 22 

E, eV 

FIG. 2: Fermi-LAT data (black circles) for the EGRB and 
UHECR data from HiRes (dots) together with UHE neutrino 
(stars) and photon (boxes) fluxes for £ max = 10 21 eV, z max = 
2, m — and a g = 2.0 (blue, open) and a g = 2.6 (red, filled 
symbols). 



10° 
10 5 
10 4 t 



o 10 
> 

to 



LU 



LU 



10' r 




BAIKAL ANITA-lite 
L' cascade U9er ^- JEM-EUSO „ a A^ 



IceCube 5yr 




10 



FIG. 3: Upper limits on the all-flavor UHE neutrino flux 
and expected sensitivities [3 together with the cascade limit 
("_E~ 2 cascade"). Also shown are realistic fluxes of cosmo- 
genic neutrinos marked by their spectral index a g — 2.6 (dip 
model) and a g = 2.0 (ankle model) together with neutrino 
fluxes optimized for detection by IceCube and JEM-EUSO 
(as described in Section IV), which marked in the Figure by 
their respective -E max values in eV (10 20 and I0 22 ). 



is shown for the ankle model with a transition from galac- 
tic to extragalactic cosmic rays at 5 x 10 18 eV for the same 
values of E max and z m ax- From Fig. [T] one can see that 
both models are allowed by the cascade limit. 

The MC simulation allows us to test the universality 
of the cascade spectrum. If a cascade is initiated by a 
photon or an electron of very high energy, the energy 
spectrum of the resulting cascade photons depends only 
weakly on the energy of the primary particle for a suffi- 
ciently large number of cascade steps. This universality is 
obviously broken for the primaries injected close enough 
to an observer, if the distance is of the order of the ab- 
sorption length (see Eq. [2} . In Fig. [T] we plot the MC 
cascade spectra with a g = 2.6 and a g — 2.0 normalizing 
them by the highest energy point of the Fermi spectrum 
(red and blue stars in Fig. [TJ. The comparison of the 
three theoretical spectra at energies below the minimal 
absorption energy e a shows that the cascade spectrum is 
indeed quite universal. The shape of the cascade pho- 
ton spectra from the Monte Carlo simulation agrees rea- 
sonably well with the one analytically calculated, with a 
somewhat harder photon flux obtained with the Monte 
Carlo method in the plateau region, J(E) oc _E~ 195 . As 
a result, the maximal cascade energy density w™f x ob- 
tained using the Monte Carlo simulation is 30% smaller 
than in the analytic calculations. 

In Fig. [2 we show the obtained UHECR, neutrino and 
photon fluxes together with data from HiRes and Fermi- 
LAT for the two cases a g = 2.0 (blue) and 2.6 (red). We 
use again £ max = 10 21 eV, z m ax = 2 and normalize the 
UHECR results to the HiRes observations. While the dip 
model fits the HiRes data with x 2 = 19-5 for d.o.f. = 19, 
the ankle model cannot explain the HiRes data below 1 x 
I0 19 eV without an additional component. Clearly, the 
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FIG. 4: Photon fluxes from the Monte Carlo simulation for 
different magnetic field strengths B — 0.01 and I nG with 
Emax = 10 21 eV, z max = 2, m = and a g = 2.0. 



dip scenario without evolution and with modest values 
of -E max and z max is well compatible with the Fermi data 
(see Fig. []}. The ankle scenario with a g = 2.0 has a 
lower flux of cascade gamma-radiation and is viable too. 



IV. THE CASCADE BOUND ON UHE 
NEUTRINOS 

This is the most general bound on the UHE neutrino 
flux, based only on the production of electromagnetic 
cascades, which inevitably accompany the_production of 
pions responsible for the neutrino flux 0, LU|. It is based 
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on the approximate equality of the total energy release 
to neutrino radiation (through p7 — > — > v and to the 



cascade radiation (through pj 



and pj 



7T — > 7) in pion production process. The upper limit 
on the integral flux J„(> E) of neutrinos of all flavors is 
given by the following chain of inequalities, 

>Z» > — / E'J v {E')dE' > — EJ U (> E), 
c Je c 

where x and w^ as are the energy density of the cascade 
radiation allowed by the Fermi data and that produced 
only by pions, respectively. For the sake of comparison 
with experimental upper bounds, where a E~ 2 neutrino 
spectrum is usually assumed, we give the upper limit for 
the differential cosmogenic neutrino flux of three neutrino 
flavors with a E~ 2 spectrum and as function of the ratio 
of energy densities of pair- and pion-produced cascades 

cas / cas ' 



only 22% of the cascade energy is lost due to synchrotron 
radiation. 

The case considered above corresponds to the E~ 2 up- 
per limit shown in Fig. [3J For cosmogenic neutrinos pro- 
duced by protons interacting with CMB the fraction of 
energy lost from the cascade is less because of the steeper 
generation spectrum Q P (E) cx E~^ g with (3 g ~ 2.3 — 2.7. 
As a result the cascade energy is produced mainly by 
low-energy electrons, for which IC dominates. The ratio 
([7]) is given now by 



1 - (E?/E min ) 



-(A, 



l-{E max /E min )-^- 2 )' 



(8) 



and is for all practical cases very close to 1. Our photon 
fluxes from MC simulations with B = 0.01 and InG 
displayed in Fig. 0] show indeed only minor differences, 
compatible with the analytical estimate given above. 



E 2 J V {E) < 



47r ho.(E max /E min ) 1 



'Ma, 



(6) V. COSMOGENIC UHE NEUTRINO FLUXES 



This limit is plotted in Fig. H as a red line labeled 
'_E~ 2 cascade' together with existing upper limits from 
various experiments and the expected sensitivity of Ice- 
Cube and JEM-EUSO 0, []J]. Equation © gives the 
general upper limit on the neutrino flux using the E~ 2 
assumption. However, each particular model for cosmo- 
genic neutrinos can be checked for consistency with the 
Fermi bound straightforwardly, as described in Section 
II. Namely, w ca s can be calculated from Eq. ([5]) and com- 
pared with w™f x = 5.8 x 10~ 7 eV/cm 3 . 

We discuss now the impact of magnetic fields on the 
cascade limit. In the presence of magnetic fields, the en- 
ergy of electrons is partly dissipated in the form of syn- 
chrotron radiation. The critical energy E™ of electrons 
above which synchrotron energy losses dominate is de- 
termined by the relation (dE/dt) syn — {dE/dt)\Q 1 where 
the indices 'syn' and 'IC are related to synchrotron and 
IC losses, respectively. 

In the case of a single electron with energy Eq > E™, 
the usual electromagnetic cascade is suppressed until the 
electron energy drops below E". In this case, the total 
cascade energy is reduced by the factor E" /Eq compared 
to the case without magnetic field. However, the cascade 
is initiated by many electrons with production spectrum 
oc E~ 2 , as required for a -E~ 2 -upper limit. Then the 
cascade energy density u> cas {E e )dE e is proportional to 
E v J v (E v )dE v , and the ratio of the cascade energy den- 
sity in presence of a magnetic field oj^ s and in its absence 
uj cas is given by 



ln(£ e "/£ e mi ") 
lii(£ e mM /£ e mil1 ) 



For the case of a strong magnetic field B ~ 1 nG one 
obtains E P CI - 2 x 10 18 eV. Taking £ roax - 1 x 10 21 eV 



and E„ 



1 x 10 9 eV, we find 



We discuss now the cosmogenic neutrino fluxes com- 
patible with the two conditions that the parent proton 
fluxes provide a good fit to the HiRes data and that the 
resulting EGRB respects the cascade bound. The latter 
is imposed by requiring that w cas calculated with the help 
of Eq. ([5]) is smaller than w™g X deduced from the Fermi 
data. 

We begin with our reference models, given by the dip 
(a g = 2.6) and ankle (a g = 2.0) model normalized to the 
HiRes data, and using E max = 1.0 x 10 21 eV, z max = 2 
and non-evolution. These are conservative models which 
give in case of a g — 2.6 — 2.7 the lowest neutrino fluxes for 
the proton-dominated mass composition. These fluxes 
labeled in Fig. [3] by the generation indices 2.0 and 2.6 
are shown in the central part of the figure. They are 
undetectable by Auger and the planned detector JEM- 
EUSO. 

Next we extend our reference models, allowing larger 
£ max and cosmological evolution. This results in higher 
neutrino fluxes, limited however still by w C as- The bounds 
on the parameters of these models are shown in Fig. [5] 
The two panels of this figure show for E max = 1 x 10 21 eV 
(left panel) and -E max = 1 x 10 22 eV (right panel), how 
the two parameters describing source evolution, m and 
z mM , are limited by w max = 5.8 x 10 -7 eV/cm 3 . Gener- 
ally, strong evolution (m, z max ) with m > 3 and z max <L 4 
is excluded, in accordance with the cases presented in 
Table 1. The evolution in the dip model (a g = 2.6) is 
restricted stronger, because of the increased contribution 
from e + e~ pair-production. 

The cosmogenic neutrino flux can become detectable 
only in the case of source evolution and large -E max . Two 
extreme models of such neutrino fluxes are shown in the 
lower-right and lower-left corners of Fig. [3] Both mod- 
els respect the bound from the observed UHECR flux 
0.78, i.e. and from w™| x derived from the Fermi-LAT data, but 



(7) 
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FIG. 5: Range of allowed evolution parameters, m and z max , for extended reference models with fixed E maK = 1 X 10 21 eV (left 
panel) and -E ma x = 1 X 10 22 eV (right panel). The cascade energy density aj cas is shown as function of m by the solid lines for 
the ankle model (a g = 2.0), and dashed lines for the dip model (a g = 2.6). The numbers on the lines show z max . The allowed 
parameters correspond to part of the curves below w^f* = 5.8 x 10 -7 eV/cm 3 shown by the red horizontal line. 



use extreme values for the model parameters. Choosing 
the parameters for the model in the lower-right corner 
(the curve marked 10 22 ) we try to reach the sensitivity 
of JEM-EUSO. Since a soft spectrum increases w C as, we 
choose the hard spectrum with a g = 2.0, while E max 
should be as large as possible. By other words we search 
for the extension of the ankle reference model with al- 
lowed evolution and large -Emax- We choose E max = 
1 x 10 22 cV, with z max = 2 and evolution parameter 
rn = 3. Normalized to the HiRes data, this model has 
^cas = 3.3 x 10~ 7 eV/cm 3 , i.e. is somewhat below the cas- 
cade limit (see also Fig. [5]). For such values, the neutrino 
flux is marginally detectable by JEM-EUSO. 

In the lower-left corner (the curve marked 10 20 ) we aim 
to cosmogenic neutrino detection by IceCube. Here we 
should increase the low-energy tail of the neutrino flux 
and suppress the pair-produced cascade radiation. To 
that end, we use a g — 2.0 with strong evolution to en- 
hance the flux of low-energy neutrinos. The maximum 
acceleration energy can be low, e.g. E max = 1 x 10 20 eV. 
Moreover, we choose evolution with m = 3.0 and z max = 
6.0, which results in cj cas = 5.5 x 10~ 7 cV/cm 3 w w™| x . 
As our calculations show, the flux is only marginally de- 
tectable by IceCube even for these extreme parameters. 

The two models above demonstrate that even for ex- 
treme assumptions cosmogenic neutrinos remain unde- 
tectable by existing detectors such as Auger, and could 
be only marginally observed by IceCube and by future 
detectors JEM-EUSO and Auger-North (with sensitivity 
to neutrinos 5-6 times higher than Auger-South). 

The observation of radio emission from neutrino- 
induced air showers provides an effective method for the 
detection of low fluxes of cosmogenic neutrinos from the 
highest energy part of their spectrum. The upper limit 
on UHE cosmogenic neutrino flux from the most restric- 
tive experiment of this type, ANITA, is shown in Fig. 3 



(Gorham et al. [HI). Recently, several particles with 
energies above 1 x 10 19 eV have been detected there 
[20l | . The high energy threshold is a disadvantage of this 
method. In the recently proposed ARIANNA detector 
[2lj . the threshold might be lowered to about 10 17 eV 
while monitoring 900 km 2 of Antarctic ice. 

A very sensitive instrument for UHE neutrino detec- 
tion has been proposed in the project LORD (Lunar Or- 
bital Radio Detector) (22|, where a detector on a lunar 
satellite can observe the neutrino-produced radio-signal 
from lunar regolith. The sensitivity of this instrument, 
as estimated by the authors of the project, should be suf- 
ficient for the measurement of the cosmogenic neutrino 
fluxes shown in Fig. [3] by curves 10 21 . 

Before concluding, we would like to compare the re- 
sults of this investigation to the ones of Ahlers et al. [23[ 
that appeared after ours in the arXiv. While the main 
goal of our work was to derive an upper limit on the cos- 
mogenic neutrino flux, the authors of Ref. [23[ aimed at 
exploring the allowed parameter space of UHECR mod- 
els, notably of those predicting maximal neutrino fluxes. 
These authors used as their criterion for the rejection 
of UHECR models w™f x = 5.8 x 10~ 7 eV/cm 3 from our 
calculations, and thus the derived maximally allowed cos- 
mogenic neutrino fluxes should coincide. The largest cos- 
mogenic neutrino fluxes presented in Fig. 4 of Ref. [23| are 
very similar to our fluxes obtained in the extreme mod- 
els with strong cosmological evolution (e.g. the curve 10 22 
in Fig. [3]), both exceeding our reference cases (a g = 2.6 
and a g = 2.0 without evolution) by an order of mag- 
nitude at E ~ 10 18 10 19 eV. It is noteworthy that a 
much stronger cosmological evolution was considered in 
the calculations of Ref. [2J|. Among other differences, the 
authors of Ref. [23| assumed that the IceCube sensitiv- 
ity extends up to 10 19 eV, while we used E max = 10 17 eV 
following Ref. [H|. 
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VI. SUMMARY 

Wc have used a recent measurement of the EGRB by 
Fermi-LAT to constrain models for UHECR and cosmo- 
genic UHE neutrinos and to demonstrate that the latter 
are not detectable with the present experimental sensi- 
tivity. Both the dip and ankle model without or with 
weak evolution are consistent with the Fermi-LAT mea- 
surement of the EGRB. The cosmogenic neutrino flux is 
strongly limited by the new upper cascade bound and 
undetectable for a conservative choice of parameters by 
Auger-North and JEM-EUSO. Only for an extreme set 
of parameters, E max > 1 x 10 22 eV and w ca s ~ w™f x , 
the cosmogenic flux is marginally detectable by JEM- 
EUSO. To achieve the observation of cosmogenic neutri- 
nos for less extreme parameters, the detection threshold 
of JEM-EUSO (in the tilted mode) must be lowered down 



to 1 x 10 19 eV and the sensitivity of Auger-North should 
be increased by factor ~ 20 in comparison with Auger- 
South. The further development of radio-detection meth- 
ods gives another hope for detection of small fluxes of 
cosmogenic neutrinos. 

The results of our paper emphasize the necessity to 
develop more sensitive methods for the detection of cos- 
mogenic neutrinos. 
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